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Weakly interacting massive particles (WIMPs) are among the favored candidates for cold dark 
matter in the universe. The phenomenology of supersymmetric WIMPs has been quite developed 
during recent years. However, there are other possibilities which have not been discussed as much. 

One example is a right-handed massive neutrino, which has recently been proposed in the context 
of a version of the Zee model for massive neutrinos. This TeV-scale, leptonic WIMP (or LIMP, 
for short) may at first sight appear to be essentially undetectable. However, we point out that 
the radiatively induced annihilation rate into leptons and photons is bound to be substantial, and 
provides a conspicuous gamma-ray signature for annihilations in the galactic halo. This gives a 
window of opportunity for Air Cerenkov Telescopes with ability to observe the galactic center, such 
as the HESS and CANGAROO arrays, and also for the GLAST space telescope. In addition, the 
contribution to the positron cosmic ray flux is in principle detectable, but this would require very 
strong local density enhancements in the dark matter halo distribution. 

PACS numbers: 95.35.+d, 14.60.St, 95.85.Pw, 95.85.Ry, 98.70.Rz 


I. INTRODUCTION 

Recently, Krauss, Nasri and Trodden ([Q, KNT in the following) proposed an interesting model, where a right- 
handed neutrino of mass on the order of a few TeV plays a crucial role in giving mass to the otherwise massless 
standard model neutrinos through a high-order loop mechanism. This is a version of the Zee model Q, which has 
been quite successful is reproducing the observed mass and mixing pattern of solar and atmospheric neutrinos (see, 
e.g., ||). The particle content of the Zee model is given by two Higgs doublets and d> 2 , and a charged field S 
which transforms as a singlet under SU( 2), with Lagrangian 

£ Zee = fapL^CiT2LpS + + /J.^iT 2 ^ 2 S~ + h.C. (1) 

KNT consider a variant where neutrino masses appear only at the three loop level. To achieve this they supplement 
the SM fields with two charged singlet scalars Si and S 2 and one right handed neutrino Nr. Lepton number is broken 
explicitly by including a Majorana mass term for the right-handed neutrino, and imposing a discrete Z 2 symmetry 
under which the SM fields and Si are singlets but S 2 and Nr transform as 

Z 2 :{S2,N r }^{-S 2 ,-Nr} , (2) 

forbidding Dirac masses for the neutrinos. This gives the Lagrangian 

^KNT = + QciNrSzIoir 

-YMrNrCNr + V(Si, S 2 ) + h.c. , (3) 

in which the potential V(Si, S 2 ) contains a (S\S%) 2 coupling. It is assumed a mild hierarchy of masses Mr < Mr 1 < 
Ms 2 ~ TeV and that the Yukawa couplings f a/ s, g a are of order unity, making Nr stable in view of the discrete 
symmetry. Left-handed Majorana neutrino masses are induced at three-loop order. For M$ 2 ~ TeV, KNT find an 
effective dimension-five effective mass scale of A > 10 9 GeV, giving neutrino masses at the 0.1 eV scale without 
involving fundamental mass scales significantly larger than a TeV. 
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II. TWO-BODY TREE-LEVEL ANNIHILATION RATES 


The discrete symmetry and the fact that Nr is lighter than the charged scalars means that Nr becomes stable and 
therefore a natural dark matter candidate. Through S 2 exchange it coupled to charged leptons in the early universe 
strongly enough to give the correct relic density, but extremely weakly today. Current direct dark matter detectors 
employ scattering on nucleons, and will not be sensitive to the available leptonic interactions. One could imagine 
techniques based on atom level transitions Q, but a simple estimate shows that the rate is typically less than one 
event per ton per year, thus several tons of active detector material would be needed - a very difficult task given that 
a gaseous phase detector will probably be necessary. 

Similarly, indirect detection through neutrinos from the Earth or the Sun will not be possible, since the cross section 
for capture is negligibly small. 

We therefore turn to indirect detection through annihilation in the galactic dark matter halo. First, it may be 
useful to review the thermal production mechanism in the early universe. Since the Nr is a Majorana particle, its 
annihilation into a lepton pair shows the usual helicity suppression av oc m 2 for the S-wave. Repeating the original 
calculation of Goldberg for photinos |Q, the cross section can be written 


av ( NrNr - £+£-) = 


at 


+ / 2 ) 2 


m t + o 


1 + / 4 
3 V(l + / 2 ) 2 


2 2 
m N v 


( 4 ) 


where me is the charged lepton mass, ms = frriN the S2 mass and to at the LIMP (Nr) mass. Note that for / not 
too much larger than unity, the factor in square brackets is typically close to 0.5. This means that for m at in the TeV 
range, the P-wave cross section (the term proportional to v 2 in Eq. (|J)) will determine the freeze-out temperature Tf 
and therefore the relic abundance |Q. As usual for a WIMP, one finds Tj/toat ~ 1/20, and so ( v 2 )/ = 6Tf/mN ~ 0.3. 
The requirement that the LIMP be the dark matter particle, with relic abundance fl^h 2 ss 0.1, as indicated from 
a joint analysis of cosmic microwave background and large scale structure data (see, e.g., Q), then fixes the cross 
section at freeze-out to be |?j, |f| 


E ( av )f 


Yte at _ 

8071771^(1 + f 2 ) 2 


3 x 10 “ 26 cm 3 s _1 . 


( 5 ) 


For a given Nr of mass to at, this fixes the normalization of the combination 4“ / 2 ) 2 which appears in 

many annihilation formulas. We can then estimate the total annihilation rate into lepton pairs £ + £“ at rest (i.e. in 
S-wave) putting v = 0 in Ecp ([|). We find 

av (NrNr —*■ £ + £ _ ) « 3 x 10 ~ 25 f 2 ^ cm 3 s _1 = 10 -25 (—§-\ cm 3 s -1 (flavor universal). (6) 

\l^e ae m N/ \ m N) 

In the remainder of this paper we will assume flavor universality, namely g e = = dr = d- For a TeV LIMP, 

this looks like a phenomenally small annihilation rate for the e + e _ and g + g~ channels. Even for t + t“, the helicity 
suppression is of the order of 10 -6 (at that level, the P-wave term in the galactic halo starts to contribute since typical 
galactic velocities correspond to v ~ 10 -3 ). One would therefore be tempted to conclude that the LIMP, despite being 
a good dark matter candidate, has little chance to be detected in any direct or indirect detection experiment. We 
note that in order that the couplings ge not be much larger than unity, for nearly degenerate Nr and S 2 , mN cannot 
be too much heavier than 1 TeV, but could be significantly lighter. 


III. THE DETECTABILITY OF Nr THROUGH GAMMA-RAYS 


The small annihilation rate, unlike the small scattering rate (which is caused by the leptonic nature of the candidate), 
is due to the Majorana nature of Nr which implies the absence of large S-wave two-body annihilation at tree level [Q. 
Any higher order process which does not have this helicity suppression will have a chance to dominate the annihilation 
rate completely. This is precisely the kind of interesting situation investigated many years ago for the case of a pure 
photino coupling through light selectrons to electrons, positrons and photons j^]. (That particular candidate is now 
less plausible in view of results from the LEP accelerator.) In particular, we can immediately take over the results of 
jP] to write the differential cross section for the radiative, non-helicity suppressed process NrNr — > as (here 

the annihilation rate T = av/ 2 ) 


dT 


(NrNr -> £+£" 7 ) = 3 X 10“ 


■25 dt a(m 2 N +m 2 s ) 2 

J2e at mJl N 


F(E 1 ,E t ) cm 3 s _1 GeV 
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( 7 ) 
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where we have neglected the lepton masses (a good approximation also for the for TeV LIMPs), and where 

( rn,N ~ Ej)(2m% - ArriNEf - 2toat£' 7 + 2 Ef + 2 EfE 1 + E 2 ) 

7 ’ 1 (3rn 2 N — 2mj<iEi — 2m^E 1 + rn 2 s ) 2 {m 2 N — 2mjsiEf — m 2 s ) 2 

Note that the total rate for gamma rays of this process is fixed by the cross section giving the relic density - it is 
independent of the values of the individual gf in this massless limit. The strength of this process is seen to be of the 
order of a/V times the annihilation rate at freeze-out, which is orders of magnitude larger than the helicity-suppressed 
two-body S-wave annihilation rate. This will mean, as we shall see, that there is a hope of detecting the LIMP in 
gamma-rays, and perhaps also in an anomalous positron component of the cosmic rays. On the other hand, since all 
three charged SM leptons are too light to decay into baryons, indirect detection through antiprotons (or antideuterons) 
is not expected for the LIMP. A nice feature of the result (Eq. |?j) is its model independence. For any dark matter 
candidate of this type (a Majorana particle coupling mainly to leptons through charged scalar exchange) we expect 
this strength of the radiative annihilation signal. This is to be contrasted with the situation in the MSSM, where the 
predicted gamma-ray signal for a SUSY WIMP of any given mass depends strongly on a large number of additional 
supersymmetric parameters. However, there are astrophysical uncertainties which make the expected absolute fluxes 
difficult to estimate despite this robustness of the particle physics properties. 

By integrating Eq. |z]over Ef, (with lower and upper integration limits to at — Ej and tojv, respectively) we get the 
differential photon spectrum. Similarly, by integrating over E 1 between mjy — Ef and rriN the differential lepton 
spectrum (and the identical anti-lepton spectrum) is obtained. Apart from the overall normalization, which is fixed 
by the relic density and contains a factor (tojv /ms ) 4 strongly favoring a scenario with at most a mild hierarchy of 
these masses, the shapes of these distributions depend only on the scaled energies E/tojv and the ratio TOat/tos. 

Since all leptons can be treated as massless, the distributions are identical for e±y, /i ± 7 and so from the point 
of view of this gamma-ray signature, it is of no importance whether the Nr couples universally to leptons or not. 
There are in principle a significant number of hard photons in the tau decay chain from pion decays, but we have 
checked that the direct photons from £ ±r y strongly dominate in all cases we illustrate. For, e.g., the positron signal 
there may be differences, since muons and taus give positrons with softer spectra in their decays, as discussed in the 
next Section. In Fig. [I] (a) and (b) we show examples of these distributions for mjv/ms = 0.8 and 0.2, plotted on 
logarithmic and linear scales, respectively. As can be seen, both the gamma and the lepton spectra are exceptionally 
hard, peaking near the maximum energy (which in both cases is equal to uin). This is optimal from the point of view 
of detection, as most conceivable sources of background tend to give soft spectra, rapidly falling with energy. 




FIG. 1: (a) The differential photon spectrum for the process NrNr —> £ + £~-y, normalized to unity, for mjv/m-s = 0.8 (solid 
line) and wat/tos = 0.2 (dash-dotted line), as well as the lepton (or the identical anti-lepton) spectrum for the same mass 
ratios (dotted line and dashed line, respectively), (b) Same as in (a) but plotted using a linear scale on both axes. 

The obvious place to search for the gamma-ray signal is in the direction of the galactic center, where the dark 
matter distribution is expected to be strongly enhanced. N-body simulations indicate that cold dark matter, of which 
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the LIMP is an example, is even likely to form a density cusp p(r) ~ 1/r near the center (the so-called Navarro-Frenk- 
White or NFW profile jpj). The observational situation concerning rotation curves and the distribution of dark 
matter in spiral galaxies, including our own, is however far from clear at the moment, and the same is true for the 
N-body calculations where different groups still seem get different results mm- Even more unclear is the interesting 
pos sibility of the massive black hole transforming the cusp into a very sharp spike with dramatically enhanced density 

iii- 

The best we can do at the moment is to assume an NFW profile, keeping in mind that the predicted rates can be 
either much lower (in the case of the absence of both a cusp and a spike) or much higher (in the case of a steeper 
cusp, as some N-body simulations indicate, or a spike, or a clumpy halo [jH| as also indicated by simulations p7[). 
Here we note the importance of the dumpiness of the halo: since the annihilation rate is proportional to the square 
of the density p, it is enhanced over the prediction for a smooth halo by a dumpiness parameter C = (p 2 )/{p) 2 > 1. 

The rate of gamma-rays from a given halo profile has been elaborated in ]li|. Combining the results of that work 
with those given here, we predict a gamma-ray flux from the galactic center (treated as a point source, so the units 
are in cm -2 s ” 1 GeV” 1 ) 


cd> 7 

dEj 


= 1.3 x 1CT 11 


/ 1 TeV \ 2 /(J) 

V J V 100 / 


( m N + m |) 2 

m N 


/ m n —E^y 


dEiF(E^, Eg), 


(9) 


where the line-of-sight integral J is defined in |h|, and (J)aq. is the average of this line-of-sight integral over the 
acceptance AH of a gamma-ray telescope; for a NFW profile, the maximal value of the product (J)ahAH ~ 100 
for AH ~ 10 ” 3 [^ 8 | (similar, or larger, values would be obtained for the Moore profile jyj p(r) ~ 1/r 1 ' 5 ). The 
background is difficult to estimate (in fact, much of the presently perceived diffuse “background” may turn out to be 

M t of the kind of signal discussed here). For the sake of the argument, we will take the simple estimate given in 
, extrapolated from the flux measured by EGRET towards the galactic center: 


d <f>bkg 

dEjdfl 


= 4.8 x 10 ” 13 


g.c. 


(I TeV 


2.7 

cm ” 2 s ” 1 GeV ” 1 sr” 1 . 


( 10 ) 


An inevitable consequence of the LIMP coupling to a charged scalar is that the loop-induced process N R N R —> 77 
should also occur |J, with cross section 


av {N R N R -> 77 ) 


'} x 4 x 10” 3 ° cm 3 s” 1 . 


(ii) 


Note that since all leptons contribute coherently to the loop amplitude, the rate is more favorable for a universal 
lepton coupling. In the case g e = = g T = g , the value for av becomes 1.2 x 10 ” 29 cm 3 s ” 1 which is quite sizable. 

Note that this result is independent of the N R and S '2 masses jf|. (One may also note that since S 2 is an SU(2) 
singlet, there is no corresponding Z 7 final state.) 

The line flux then becomes 


$® 7 '(.E 7 = tojv) ~ 1.2 x 10 13 T e V \ cm 2 g 1 (12) 

” V m N J 

This would appear as a sharp spike (with relative width of the order of the Doppler broadening due to LIMP motion, 
i.e. of the order of 10” 3 ), in an instrument with perfect energy resolution. However, ACTs have an energy resolution 
of the order of 5 — 10% at best. If we assume 5%, the gamma spike for a TeV LIMP spreads out over 50 GeV, a 
differential rate near a TeV similar to the radiative processes discussed above. 

For smaller LIMP masses of the order of 100 GeV, one gets into the mass range where the GLAST gamma- 
ray telescope @ will be operative (it will detect gamma-rays up to 300 GeV). The energy resolution for GLAST 
will depend on the angle of incidence and may reach a few percent, making it possible for the line to stand out 
against background. The detection of such a line would eliminate all possible confusion with plausible astrophysical 
backgrounds. 

In Fig. [| (a) we show the flux predicted from Eq. (|j) and Eq. (|5 a]) together with the background estimate Eq. ([it]) 
for a 100 GeV LIMP, and 110 GeV S 2 , NFW profile and AH = 10” 3 . For this energy range, we have used an energy 
resolution of 3%. It may be difficult to push the N R mass much below 100 GeV without fine tuning the parameters of 
the KNT model (and the S 2 mass is also bounded by LEP results to be larger than around 100 GeV). We note that 
with the S 2 (and thus Si) only 10% higher in mass than the N R , a careful analysis should really take coannihilations 
with S 2 into account when estimating the relic density, and thereby fixing the interaction strength. In analogy to 
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SUSY models where R-parity stabilizes the lightest superpartner, in this model the Z 2 symmetry stabilizes the lightest 
odd state; a complete calculation would include all odd states, in this case N R and S 2 . However, coannihilations are 
unlikely to change the relic density by more than a factor of a few, a small correction compared with the very large 
astrophysical uncertainties related to the radial distribution of dark matter near the galactic center. 

The natural mass range for the LIMP is around 1 TeV, where GLAST runs out of sensitivity but where ground- 
based arrays of Air Cerenkov Telescopes with large collecting area can detect a signal. Indeed, there are already 
such arrays of telescopes planned or in operation such as CANGAROO [ pT| , HESS fffi] , VERITAS [^3) and MAGIC 
fj4| . In particular, CANGAROO and HESS are well located to observe the galactic center for a sizable fraction of 
their observing time. As can be seen from the figure, the signal with these assumptions would stand out from the 
gamma-ray background. (We do not enter here into the more technical issue of rejecting other types of background, 
such as from hadrons and electrons, where there is a steady improvement in the techniques employed.) 





FIG. 2: (a) The total gamma-ray flux expected from a AS2 = 10~ 3 sr cone around the galactic center (solid line). The flux is 
composed by a power-law background extrapolated from EGRET data (dotted line) and a 100 GeV LIMP annihilating with 
a cusped (NFW) density profile through a 110 GeV scalar S 2 , giving both a continuous spectrum and a 2y line. An energy 
resolution of 3% has been assumed for the line signal, (b) Same as (a) for a 1 TeV LIMP, ms 2 = 1-1 TeV. Here the line has 
been smeared by an assumed energy resolution of 5%. (c) Same as (b) for an 8 TeV LIMP, ms 2 = 8.8 TeV. 

In Fig. H (b) and (c), the curves are shown for a LIMP of mass mjv = 1 TeV, ms 2 = 1-1 TeV, and mw = 8 TeV, 
171 s 2 = 8.8 TeV, respectively. These should be clearly observable with a very conspicuous “bump” in the spectrum, 
for the halo parameters chosen. We note with interest that preliminary results from the CANGAROO collaboration 
indeed show an excess flux of TeV gamma-rays from the galactic center |2l| . The absolute flux level for this possible 
signal seems higher than that predicted in Eq. (!), so an enhancement beyond that provided by the NFW profile 
would then be indicated. 

An interesting question is what would happen to the leptons. To be specific, assume that there is lepton universality 
in the Nr couplings. Then there will be equal amounts (apart from very small lepton mass corrections) of ^ and 
er- leptons. The produced and will decay into lower-energy electrons and positrons (see next Section), and the 
primary high-energy electrons and positrons will quickly radiate due to synchrotron and inverse Compton processes. 
In fact, an analysis by Bertone, Sigl and Silk J2lj] points to the existence of a radio signal that may be caused by 
just the kind of TeV-scale dark matter particles discussed here, if the magnetic field near the galactic center is strong 
enough. In any case, this would be a signal to search for, were a TeV gamma-ray excess to be confirmed. 


IV. PREDICTIONS FOR COSMIC-RAY POSITRONS 

In the previous section the annihilation process N R N R —* l + l ~7 was discussed with the idea of detecting the high 
energy photons produced. We can also make a prediction for the cosmic ray positron flux from this process. As shown 
below, the positron flux is quite low, and not likely to be detectable unless the dumpiness factor of the galactic halo 
is quite large. 

Positrons might come from any of the three leptonic channels. In the direct process N R N R —1 e + e _ 7 the positrons 
are quite energetic. The energy is downgraded somewhat in the muon channel, since some of the energy goes to 
neutrinos: N R N R -4 /i+/x“ 7 —> e + e~v e v e v^. The spectrum is softest for the tau channel N R N R —> t + t~ 7 , since 








6 


there are many hadronic modes with longer decay chains yielding positrons (and electrons) in the end. We also include 
the simpler processes N R N R —> i + t~, For electrons this is negligible unless M N is very small and Mg 3> M N is much 
larger. However, for muons and especially taus, the simple process can be important as the helicity suppression is 
less severe. All of this can be accounted for with standard tools; in particular we have used results from the PYTHIA 
event generator ^7| as tabulated in DarkSUSY Q for the positron spectra from the muon and tau decay chains. 
Folding in the annihilation cross section, we can then compute the volume production rate of positrons in the Galactic 
halo ( dT/dE ) in units of cm 3 s _1 GeV -1 . 

The propagation of positrons in the Galactic environment is complex. The gyroradii of charged particles in the 
tangled magnetic fields are quite small, and the motions can be modeled as a diffusion process. Furthermore, positrons 
lose energy rapidly to both synchrotron radiation, and in addition to inverse Compton scattering with the cosmic 
microwave background and with diffuse starlight. We use the propagation model in Ref. to calculate the observed 
flux at the Earth. The model includes diffusion in an infinite slab, a reasonable approximation given the fact that 
detectable positrons are produced within a few kpc, thus an outer radial boundary is unimportant. The energy losses 
to synchrotron and inverse Compton processes are included. This model is in rough agreement with earlier work [ po[ , 
though the inclusion of inverse Compton scattering from starlight doubles the energy loss rate. More sophisticated 
models give similar results |n] . We note here that the question of the halo profile is much less important in calculating 
the positron flux, as their effective range is only a few kpc. Thus, it is the less uncertain local dark matter density 
that is important for the positron flux. 

The HEAT collaboration has in fact reported an excess of positrons above about 10 GeV in the cosmic ray positron 
fraction / = e + /(e + + e ~) Q, later confirmed in a different instrument by HEAT-pbar (33). We now discuss the 
possibility that this excess might be due to annihilating LIMPs. The fluxes of positrons are typically too small by a 
factor of 10 4 if a smooth halo is assumed. This might seem hopeless, but the dumpiness parameter could be this large, 
boosting the signal accordingly. In order to produce an excess at 10 GeV, lighter LIMPs are preferred. Assuming that 
the LIMP must be heavier than 100 GeV to escape discovery (actually, it is the heavier charged scalars Si and S 2 
which would be seen at e.g. LEP), we illustrate some LIMP models that give an acceptable positron flux to explain 
the HEAT results in Fig. (albeit with the very large dumpiness parameter). Other explanations for this excess have 
been put forward elsewhere |29| 34|, |35|. In the second panel of Fig. [ 3 ] we show some LIMP models with larger masses 
that do not do well in explaining the current excess. However, future experiments such as AMS and PAMELA 
p8[ may be able to probe such models with higher accuracy, and we estimate that boost factors a factor of five or 
ten smaller (of order 10 3 ) would still produce an interesting signal for such experiments. 
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FIG. 3: Cosmic ray positron fraction for LIMP models. The data from the two HEAT experiments in illustrated, along with 
the expected value (dotted line), (a) Variation of Ms with values of (top to bottom) 120, 150, and 200 GeV, with Mjv fixed at 
100 GeV. The boost factor is fixed at 2 x 10 4 . (b) Variation in Mjv, with values of 100, 150, 250, 400, 600, and 1000 GeV. In 
all cases Ms = 1.2Mjv. Here a boost factor of 4 x 10 3 is used. 
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V. DISCUSSION AND CONCLUSIONS 

We have discussed the possibility of astrophysically detecting a right-handed neutrino cold dark matter candidate 
arising in a recently proposed modified Zee model. As the interactions are only leptonic, current elastic scattering 
experiments are not sensitive to this particle. Annihilations in the galactic halo, with both two and three body final 
states, may provide detectable numbers of high energy photons or positrons, though these predictions rely heavily on 
the structure of the Galactic halo, both in core profile and in spectrum of substructure. 
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